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We suggest that clusters or domains of topological charge and action density occur in the QCD 
vacuum as an effect of singularities in gauge fields and can simultaneously lead to confinement 
and chiral symmetry breaking. The string constant, condensates and topological susceptibility 
are estimated within a simplified model of hyperspherical domains with interiors of constant field 
strength with reasonable values obtained. Propagators of dynamical quarks and gluons have compact 
support in configuration space, thus having entire Fourier transforms, which gives rise to their 
confinement. 
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Since Gribov and Singer , evidence has mounted that 
singularities in gauge vector potentials are intrinsically 
unavoidable in nonabelian theories. Several scenarios of 
confinement based on "condensation" of singular confi- 
rations such as monopoles and vortices, typical for the 
maximal abclian, Laplacian or axial gauges, are the sub- 
ject of intensive analytical and numerical lattice studies 
today . These focus mainly on the singular fields them- 
selves, their emergence from gauge fixing, their topolog- 
ical properties and role in the QCD functional integral. 

Here we emphasise a less studied effect due to singular 
gauge fields, namely their restrictive influence on fluctu- 
ations in the vicinity of singularities ||,^ . The subtleties 
of separating fields into regular and singular parts and the 
behaviour of regular fields at the singularities are irrele- 
vant if one could calculate the QCD functional integral 
"exactly" . But these issues become crucial if one under- 
takes approximations In gauge invariant quantities 
singularities due to ambiguities in gauge fixing should not 
occur. This can happen in the action either due to can- 
cellations between derivative and commutor parts in the 
field strength if the singularity is topologically nontrivial 
(monopole or vortex) or due to finitness of both terms 
separately for topologically trivial singularities (domain 
walls, which, in general, can also appear j^). The can- 
cellation of singularities in the action density can be pro- 
hibited by fluctuations around them. Thus finiteness of 
the action implies specific constraints on fluctuations as 
is discussed in Ref. where Polyakov gauge monopoles 
are considered. 

We formulate a model partition function which incor- 
porates singularities in gauge fields effectively via their 
restrictive effect on fluctuations. We assume that singu- 
larities are present in general in gauge potentials, and in 
their vicinity one can divide an arbitrary fleld A into sin- 
gular pure gauge S and regular Q parts: = S*^ + Q^- 



The field strength for pure gauge S vanishes, but the 
field strength and the action density for A is singu- 
lar unless the field Q satisfies certain conditions in the 
vicinity of the singularities in S. To be explicit, at the 
cost of generality, we further assume that singularities in 
vector potentials are concentrated on hypersurfaces dVj 
{j — 1, . . . , N) in Euclidean space of volume V, in the 
vicinity of which gauge fields can be divided as above 

into a sum of a singular pure gauge Sp! and regular fluc- 

(i) 

tuation part Qji , with a colour vector n'j associated with 
gU) , Por such fields to have finite action the fluctuations 
charged with respect to rij must obey specific conditions 
on dVj. The interiors of these regions thus constitute 
"domains" Vj . Demanding finiteness of the classical ac- 
tion density, one arrives at 



for X £ dVj, with the adjoint matrix hj = T°'n°: in the 
condition for gluons, and a bag- like boundary condition 
for quarks, Tjjj^ix) being a unit vector normal to dVj. 

Equations (^ indicate that gauge modes neutral with 
respect to are not restricted and provide for interac- 
tions between domains. In a given domain Vj the effect of 
fluctuations in the rest of the system is manifested by an 
external gauge field Bj^ neutral with respect to n° . This 
enables an approximation in which domains are treated 
as decoupled but, simultaneously, a compensating mean 
field is introduced in their interiors. The model becomes 
analytically tractable if we consider spherical domains 
with fixed radius R and approximate the mean field in 
Vj by a covariantly constant (anti-)self-dual configuration 
with the field strength 

= cos^, -f t« sin^, , e {(2fc + l)V6}Lo, (2) 



where the parameter B = const is the same for all do- 
mains and the constant matrix n°t° belongs to the Car- 
tan subalgebra. Note that there is no source for this 
field on the boundary and therefore it should be treated 
as strictly homogeneous in all further calculations. The 
homogeneity itself appears as an approximation. 

A self-consistent mean field approach requires calcu- 
lation of the effective action as a functional of the mean 
field and characteristic functions of the domains. Its min- 
ima would give information about mean field character, 
shape and typical domain size. Needless to say that this 
problem has yet to be formulated. Nonetheless in Eqs. (||) 
we have already assumed the effective action favouring 
nonzero mean field strength parameter B and finite typ- 
ical size R. With this and arbitrary constant mean field, 
it can be shown that the effective action for a domain 
exhibits twelve degenerate discrete minima correspond- 
ing to (anti-)self-dual configurations and six values (for 
SU{3)) of the angle ^ associated with the Weyl group. 
There is also a degeneracy in the orientation of the chro- 
momagnetic field. The value — 7r/6 is specific for an 
ansatz with the effective action polynomial in TiB'^, but 
the period 7r/3 is universal. Since the volume of the do- 
main is finite the degenerate minima do not correspond 
to thermodynamical phases and have to be summed in 
the partition function. The partition function for the 
model is defined as 

Z ^jVlhnYlJ ^ J da, j T^Q" j V^P.Vi^, x 

where the thermodynamic limit assumes V, N ^ oo with 
the density w"^ = N/V taken finite. The fields Q^*', 
-01 and are subject to boundary conditions Eq.(0), in 
which the original singularities are effectively encoded. 
Interaction between the original domains is substituted 
by the mean field. A background gauge condition is im- 
posed. The integration measure ddi is 
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da,, ... = 



1 
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dai / dipi / d9i sin Oi x 



(3) 
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fe=0,l Q /=0 

Here ipi and 9i are spherical angles of the chromomag- 
netic field, tOi is the angle between the chromomagnetic 
and chromoelectric fields, is the angle in the colour 
matrix n^, is the chiral angle and Zi is the centre of 
the domain Vi. 

This partition function describes a statistical system of 
density composed of noninteracting clusters, each of 
which is characterised by a set of internal parameters and 



internal dynamics represented by the fluctuation fields. 
Correlation functions can be calculated taking the mean 
field into account explicitly and decomposing over the 
fiuctuations. First of all we consider vacuum character- 
istics of the system in the zeroth order of this expansion. 

The connected point correlator of the mean field 
strength tensor is found to be 

. . . B^:,„(a:„)) = . . . ,x„), 

N 

b;a^) = E ^^'^"Bi'M^ -{x- z,f/R^), 
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by explicit calculation using the measure, Eq. (|^). 
Translation-invariant functions S„(xi, . . . , Xn) are equal 
to the volume of the overlap region for n hyperspheres 
of radius R and centres (xi, . . . , Xn), normalised to the 
volume of a single hypersphere v = t:^R'^ /2. They are 
continuous and vanish if jxi — Xj\ > 2R: correlations in 
the background field have finite range 2R and have no 
particle interpretation. The statistical ensemble of back- 
ground fields is non-Gaussian since all connected correla- 
tors are independent and cannot be reduced to the two- 
point correlators. The simplest application of this is a 
gluon condensate 

A significant vacuum parameter for the resolution of 
the ?7yi(l) problem is the topological susceptibility [0. 
First we consider the topological charge density for colour 
SU{3) in this approximation 
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B^ 



where uij G {0, tt} depends on the duality of the j-th 
domain. The topological charge is additive 



Q = / d^xQix) = qiN+ - iV_), -Nq <Q<Nq, 
Jv 

where q — B^R'^ /16 is a 'unit' topological charge, namely 
the absolute value of the topological charge of a single 
domain, and N-^ (N-) is the number of domains with 
(anti-)self-dual field, N = N+ + N^. The probabifity 
distribution for topological charge reads 



VNiQ) 



Nl 



AfN{Q) ^ 

Un 2^ {N/2-Q/2q)\{N/2 + Q/2q)V' 



where Mn{Q) is the number of configurations with a 
given charge and Mn is the total number of configura- 
tions. The distribution is symmetric about (5 = 0, where 



it has a maximum for N even. For N odd the maxi- 
mum is at Q = ±q. Averaged topological charge is zero. 
The topological susceptibility is then determined by the 
two-point correlator of topological charge density: 



1287r2' 



Static quark confinement is examined via calculation 
of the Wilson loop 



WiL) = hm n / ^ / ^Tre 



Is 



Path ordering here is unnecessary as the matrices n*^ be- 
long to the Cartan subalgebra. For computational conve- 
nience we consider a circular contour in the {x3,X4) plane 
of radius L with centre at the origin. Calculation of the 
colour trace, integrating over spatial orientations of the 
mean field and positions of the domains and then taking 
the thermodynamic limit {N — > oo, v = V/N — 
gives an area law for a large Wilson loop L ^ R 



W{L) = e-^^i'+oCi), a = BfinBR^ 
f 

\ 



„ V3 f dx . 2\/3 f dx . 

3 / — sm X / — sm x 

2z J X z J X 



The function / is positive for z > and has a maximum 
for z = 1.557r. We choose this maximum to estimate the 
model parameters by fitting the string constant to the 
lattice result, 



(4) 



B = 947MeV, R'^ = 760MeV, 



with unit charge q = 0.15, density v ^ — 42.3fm ^ and 
the "observable" gluonic parameters of the vacuum 



= 420 MeV, x = (197 MeV)^ 
{a,/7r){F^) = 0.081 GeV'^. 



(5) 



This result indicates high density and strong background 
fields in the system. There is no separation of scales, 
\/BR « 1. Neither large domains nor stochasticity of 
background fields are seen here which a posteriori justi- 
fies the mean field averaging prescription in the partition 
function. The high density ensures area law dominance 
already at distances 2L w 1.5 — 2 fm. Taylor-expanding 
the integrand in the Wilson loop integral would reveal 
all n-point correlation functions of the background field, 
and arguments about the importance of a fast decay of 
correlators for static confinement |^ would then be seen 
to apply here. 

Due to averaging over self- and anti-self-dual configura- 
tions and angles there is no explicit violation of parity 
and chiral symmetry in the partition function. A spon- 
taneous breaking can be tested by estimating the quark 



condensate. To first order in fluctuations this requires 
calculation of the quark propagator with the condition 

z fi{x)e''''''S{x, y) - -S{x, y), [x ~ zf = R\ 

An analogous condition holds for {y — zY = R^ . Sub- 
stitution (with the upper (lower) sign for (anti-)self-dual 
domains) 

{ip + m)[P±Ho + P^O+Hi+P^O-H-i], 
0± - [l±n^B/\n\B]/2, P± = [l±75]/2, 

leads to equations for the scalar functions TC^: 

i-D^ + + 2CB)nc{x, y) = 5[x, y), 

where B = \n\B. Solutions vanishing at infinity would 
give a Green's function 7i_i divergent in the massless 
limit due to zero modes of the Dirac operator. The bag- 
like boundary conditions remove zero eigenvalues from 
the spectrum, and the massless limit is regular. 

To render our calculation transparent we choose y — 
z = and calculate the quark condensate at the domain 
centre. Then the general solutions for scalar Green's 
functions take the form {^(^ = m?' /2B + Q) 



He = A{x^\fi^) + C(;e- 



-Bx^/4 



M{l + ^i^,2,Bx^/2), 



where A(a;2|/Lt) is the vanishing at infinity scalar propaga- 
tor for mass 2_B/x, and the second term is a homogeneous 
solution regular aX x^ — expressed via the confluent 
hypergeometric function. The constants can be fit 
to implement the boundary condition. The terms toTYq 
and TTiHi vanish in the massless limit and do not con- 
tribute to the condensate. The nontrivial contribution 
comes from the term rriH-i. The bag-like conditions 
imply that on the boundary TL-i satisfies a mixed condi- 
tion with /' = df /d\x\ and the sign (— )+ corresponding 
to (anti-)self-dual domain, 

2eT^"TO7i_i = -2n'_^ - BR^^H-i, 

which leads to the relations 

r7i7i_i ~ — 



lim m7i„i - -^^^(Bi^V2)e-^"'/^ 



Tr5(0, 0) 
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Y.F{BRV2), 



\f 



^2 r ^^g2t--z(cotht-l)/2 

F(.).e-.-l + -y^-^3^(cotht-l) 



Note that the term in the propagator with nonzero trace 
has definite chirality correlated with the duality of do- 
main and is proportional to the zero mode of the Dirac 
operator: pP^O^ exp(-B2:V4) = 0. 



Averaging this result over a and (anti-)self-dual con- 
figurations and taking into account the a-dependence of 
the quark determinant |^] 

detS*"^ oc exp{ia / dxQ{x)} ~ cxp{±iqa} , 

J V 

with q being the unit topological charge, we get a finite 
result 



27r2i?3(i + g) 



^F{BR^/2) = ~-(228MeV)^ 



|fi| 



for B and R as determined above, Eq. (Q). Obviously the 
condensate vanishes for zero topological charge q. Final 
conclusions require complete calculation of the conden- 
sate and detailed analysis of the eigenvalue problem. 

Confinement of dynamical quarks and gluons can be 
studied in this context via the analytical properties of 
their propagators. As has been partly demonstrated for 
quarks, the propagators of quark and gluon fluctuations 
can be analytically calculated by reduction to the scalar 
problem, essentially that of a four-dimensional harmonic 
oscillator with total angular momentum coupled to the 
external field. The general solution is given by decom- 
position over hyperspherical harmonics. Qualitatively, 
the boundary conditions mean x— space propagators of 
charged fields are defined in regions of finite support 
where they have standard ultraviolet singularities. So 
their Fourier transforms are entire functions in the com- 



plex momentum plane 
of dynamical fields JlO|,p 



This manifests confinement 
A known consequence of en- 
tire propagators is a Regge spectrum of relativistic bound 
states 

The model displays the most essential features of QCD 
and, being analytically tractable, provides a framework 
for phenomenological applications. The first test in this 
direction should be a calculation of the spectrum of 
bound states via the Bethe-Salpeter framework, or within 
a hadronisation scheme analogous to The model 

preserves the structure of a relativistic quantum field the- 
ory, albeit a nonlocal one Certainly one has to be 
convinced that the qualitative arguments underlying the 
model are formally legitimate. Associating a colour di- 
rection with the boundary is a restrictive assumption. 
Somewhat akin to "abelian projection", it singles out a 
class of abelian mean fields, and "neutral" gluons are not 
confined. Incorporating topologically nontrivial singular- 
ities can improve this aspect of the model. 

Clustering of action and topological charge density 
with dominance by (anti-)self-dual fields in the QCD vac- 
uum is evident in lattice cooling algorithms, which, by de- 
sign, iterate towards (anti-)instantons [ p^ . Alternately, 
Horvath et al. [|4j suggested studying duality of clus- 
ters via the chirality of fermionic eigenmodes. Used with 
overlap fermions this shows strong correlation between 
clusters of topological charge and locations where several 



low eigenvalue fermionic modes are largest; the modes 
are also chiral This is regarded as evidence for in- 
stantons. Such correlation between chirality of fermionc 
modes and duality of domains is also evident in the model 
discussed above. Our estimates suggest that formation 
of clusters, predominantly (anti-)self-dual and with av- 
erage size 2R « 0.5fm, can have purely quantum origin 
whose explanation could require reference to the exis- 
tence of obstructions in gauge fixing rather than to the 
quasi-classical limit. 
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